DTP/98/50 
July 1998 



Saturation Effects in Deep Inelastic Scattering 
at low and its Implications on Diffraction 



K. Golec-BiernatQ and M. Wiisthoff 
Department of Physics, University of Durham, Durham DHl 3LE, UK 



(February 1, 2008) 



Abstract 

We present a model based on the concept of saturation for small and small 
X. With only three parameters we achieve a good description of all Deep Inelastic 
Scattering data below x = 0.01. This includes a consistent treatment of charm and 
a successful extrapolation into the photoproduction regime. The same model leads 
to a roughly constant ratio of diffractive and inclusive cross section. 



^On leave from H. Niewodniczanski Institute of Nuclear Physics, Department of Theoretical Physics, 
ul. Radzikowskiego 152, Krakow, Poland. 



1 Introduction 



The basic concept of saturation in Deep Inelastic Scattering (DIS) is related to the tran- 
sition from high to low as one observes in the total 7*p-cross section. This type of 
saturation occurs when the photon wavelength 1/Q reaches the size of the proton. We 
will include another aspect of saturation in our paper which is inherent to DIS at small x 
(small-x saturation). In this regime the partons in the proton form a dense system with 
mutual interaction and recombination which also leads to the saturation of the total cross 
section Both aspects of saturation are closely linked to confinement and unitarity. 
While the latter might be approached perturbatively the first is genuinely nonpertur- 
bative. The approach we choose here can be called QCD-inspired phenomenology and 
follows the line of refs.0, §,0. It is in spirit most similar to the ideas of the analysis 0. 

The basis for our approach is the fact that the photon splits up into a quark-antiquark 
pair (dipole), far upstream the proton target, which then scatters on the proton. In the 
pure perturbative regime the reaction is mediated by single-gluon exchange which changes 
into multi-gluon exchange when the saturation region is approached. The latter process 
can be interpreted as the interaction with a 'semiclassical field' [0. Most important 
for us is the fact that the mechanism leading to the dissociation of the photon and the 
subsequent scattering can be factorized and written in terms of a photon wave function 
convoluted with a quark-antiquark cross section a . In our analysis this cross section has 
the simple form a = ao{l — expl—r"^ / [ARKx))]} where r denotes the separation between 
the quark and antiquark and Rq is the x-dependent saturation scale: Rq{x) = (x/xq)^. 
The functional form of a can be different, important is, however, that for small r it is 
proportional to (colour transparency) while for large r the cross section approaches 
a constant value. The latter behaviour ensures saturation. The crucial element in our 
analysis is the assumption that the saturation scale Ro depends on x in such a way that 
with decreasing x-Bjorken one has to go to smaller distances (higher Q"^) in order to resolve 
the dense parton structure of the proton. The boundary in the (x, Q^)-plane along which 
saturation sets in is described by the "critical line", = 1/Rl{x). 

The recent low-Q^ data from HERA p, |10[ have triggered significant theoretical ac- 



tivity. In one class of models |TT, the description of the data is achieved by combining 
the nonperturbative vector meson contribution (VDM-model) with contributions based 
on perturbative QCD. Both models do not use the concept of small-x saturation. 



Another class of models |jT3|, |I4|, |l^] is based on Regge theory. The total cross sec- 



tion a'^ P is the sum of contributions from different Reggeons. In particular in ref. 
the leading behaviour is given by the sum of two Pomeron contributions with "hard" and 
"soft" intercepts. Again the concept of the small-x saturation is not implemented in these 
analyses. The claim in ref. |jl3[ is that saturation at present energies is not needed. How- 
ever, in the true high energy asymptotics something has to happen. The hard Pomeron 
as proposed in ref.|T^ would eventually overtake and strongly violate unitarity. Other 



approaches like those in refs. [|T6], have imposed a logarithmic behaviour in x from the 
very beginning and thus do not violate unitarity. 

The strategy we adopt in our analysis is the following. We determine the three free 
parameters of our model mentioned earlier, ctq, A and xq, by fits to all existing DIS data 
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for X < 0.01. We then study the obtained parameterization in the photoproduction 
region, where a non-zero quark mass is required to achieve a finite cross section. The 
quark mass is chosen such that the photoproduction cross section is in rough agreement 
with the data without having them included in our fit. We found that the effective slope 
Xeff of the cross section {a ~ interpolates between the "soft" Pomeron value 

^ 0.08 and the "hard" value ^ 0.29. We would like to point out that the "soft" value 
is simply a result of the extrapolation of our fits to the photoproduction region. After 
a first fit with only light quark fiavours we perform a second fit which includes charm. 
We found again a good description of all inclusive DIS data and in addition the correct 
relative charm contribution. The model we use can be straightforwardly applied to DIS 
diffractive processes leading to the interesting result that the ratio of the diffractive and 
inclusive cross section is roughly constant as a function of x and Q^. 

The content of the paper is the following. In section 2 we introduce the theoretical 
details of the model and discuss its qualitative features. We then derive a simplified 
parameterization which illustrates the concept of the critical line. This parameterization 
can also serve for a quick fit to the data. In Section 3 we describe our two fits and discuss 
their physical implications. Section 4 is devoted to the study of diffraction based on our 
model and section 5 contains conclusions. 



2 Models for Saturation 

For small values of the Bjorken variable x the photon wave function formalism has been 
established as useful tool for calculating deep inelastic and related diffractive cross sec- 
tions for 7*p scattering ^, |^. It allows to separate between the wave function of the 
photon which describes the dissociation of the photon into a quark-antiquark pair and 
the interaction of the quark-antiquark pair with the target. The photon wave function 
constitutes the calculable part of the process whereas the remainder is substantially in- 
fluenced by nonperturbative contributions and needs to be modelled. The corresponding 
diagram is shown in Fig. 1. We work in a frame where the photon with momentum q and 
the proton with momentum p are collinear. Accordingly, the distribution of the quark- 
antiquark pair is given in terms of z and (1 — z), the momentum fraction with respect 
to q, and the relative transverse separation r. For transverse (T) and longitudinally (L) 
polarized photons the 7*p-cross section takes the form Q 

(Xt,l(x,Q') = Jdh£dz\<ilT,Liz,r)\^aix,r^) , (1) 

where = {p + q)"^, = —Q^ and x = Q'^/(W'^ + Q^). The squared photon wave 
function ^t,l is given by 

I^T = ^ E 4 { + (1 - 'f\^' Kl (er) + m) KI (er) } (2) 

and 
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Figure 1: Diagrammatical representation of the basic process as discussed in the text. 



for the transverse and longitudinal photons, respectively. In the above formulae 

= z{l-z)Q^ + m}. (4) 

Kq and Ki are Mc Donald functions and the summation is performed over the quark 
flavours. The 7*p cross sections are related to the structure function F2 in the following 

way 

F2 (x, = Ft {x, Q2) + Fl (x, (5) 

and 

Ft,l{x,Q') = -^arA^^Q'). (6) 

The interaction of the qq pair with the proton is described by the dipole cross section 
(T(x,r^) which is modelled in our analysis. The most crucial element is the adoption of 
the x— dependent radius 

1 / \ ^''^ 

i?o(x) = — - , (7) 

Qo \XoJ 

which scales the quark- antiquark separation r in the dipole cross section 

a(x,r2) = aogif'') , (8) 

with 

r 

' = 2R^) ■ 

Qo = 1 GeV in (0) sets the dimension. The function g in (H) is not completely constrained. 
Important, however, is the quadratic rise at small f and the flattening off at large f. The 
latter behaviour provides saturation of the cross section (|l|), i.e. a^'^ = ax + cfl ^ 
const for small Q^. At small f, on the other hand, we have a simple scaling behaviour 
(colour transparency), a''' ^ ~ combined with a power-like dependence of (1) on x 

as typically observed in deep inelastic scattering. We choose the following simple Ansatz 
for the function g 

g (r^) = 1 - e-'' . (10) 

This Ansatz reminds of eikonalization. It should be mentioned, however, that a complete 
eikonal treatment requires the incorporation of a target profile function. The form (p!OD 
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would mean in this context that the gluonic density in the proton is evenly distributed over 
a certain area within a sharp boundary and zero beyond. A more sophisticated treatment 
including a realistic profile function can be found in ref.[^, ^. The corresponding result 
in these references can be roughly approximated by 

g{r^) = ln(l + f2) (11) 

and shows a logarithmic growth at large distances. 

For small z one can as well think of a logarithmic modification as is motivated by the 
single gluon exchange (see ref.[0) 

g{r^) = In (^1 + 1) . (12) 



Both alternatives (|TT]) and (|T^) can be easily implemented in our formalism which we 
present in the forthcoming sections. It turns out, however, that they are disfavoured in 
our analysis and we will not discuss them in detail. This does not mean that we contradict 
their physical implications. It might be that a different approach could be reconciled with 
the models (0) or (1121). 



In our analysis we fit the three parameters, ctq, Xq and A, of the dipole cross section (|D 
and ([T0|). Before going into the details of the fits it is illuminating to perform a qualitative 
analysis of the behaviour of the cross section (0) in our model. 



2.1 Qualitative Analysis 

In the following discussion we focus on the cross section in (|lD which dominates over 
(Jl. We also neglect for simplicity the quark mass. The important point in the qualitative 
analysis is the behaviour of Ki{er) in (^ for small values of er: 

Kiier) ~ — . (13) 
er 

For large values of e r the function Ki is exponentially suppressed. Thus in order to obtain 
the dominant contribution we perform the integration in (|^) for er < 1. 

The photon virtuality introduces the scale 1/Q for the transverse dimension of the qq 
pair. A pair is considered "small" when the condition r < 1/Q is fulfilled and "large" 
when r > 1/Q. Let us analyse the contribution to ([l|) coming from small pairs for which 



the condition e r = ^J z{l — z) r Q < lis satisfied for all values of z. In the case 1/Q <^ Rq, 
shown in Fig. 2a, the size of the small qq pairs is much smaller than the saturation radius 
and cr{r) ~ aof"^ /R^. The cross section (|I|) exhibits the following behaviour 



/ 1 \ 1 
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Figure 2: The profile of tlie dipole cross section for different Q. Tlie small arrows below 
the figure show how the indicated parameters change when Q decreases (for W'^ fixed). 



where the integration over z could be factored after the cancellation of the e factors. Thus 
for constant x the cross section (1|) exhibits the familiar short distance scaling behaviour, 



i.e the corresponding structure function F2{x,Q'^) is roughly constant in Q^. 

Let us now analyse the situation shown schematically in Fig. 2b in which the size of 
the qq pair is bigger than the saturation radius 1/Q > Rq. For the small pair contribution 
to (Tt we now obtain 



1 \ r 



•0 / . \ 2 




- /cir^l^-Jao^ + J d^\-2)-o - ^o + aohg^-^J. (15) 

Rl 

With respect to the power behaviour in the cross section can be viewed as being 
constant. The potential divergence due to the logarithm will be regulated by the quark 
mass in our full analysis. The actual value of the mass plays an important role in the 
description of the photoproduction region. 

A similar analysis can be performed for "large" pairs for which r > 1/Q. The inte- 
gration condition er < 1 is now satisfied when z < l/(r^(5^) and the z integration in (|1]) 
can no longer be factored out. It has to be done before integrating over r. The result in 
the end is the same as for small pairs. If the characteristic size of the qq dipole is less 
than Rq the scaling behaviour is obtained. For 1/Q > Rq the cross section ar is constant 
in Q^. A more detailed analysis gives logarithmic modifications. 

So far in our discussion we have assumed a constant saturation radius which allows a 
smooth transition of ctt between the scaling region of large Q^ and the saturation region 
of low Q^ (low-Q^ saturation). The main feature of our model, however, is the fact that 
the saturation radius Rq depends on x {Ro{x) ~ x^/^ with A > 0). In this way we have 
introduced another kind of saturation which can be called the small-x saturation. In 
terms of the parton picture it is closely related to the saturation of the gluon density |19 



An important consequence is that for fixed W the radius Rq becomes Q dependent which 
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makes the saturation more dynamical. As illustrated in Fig. 2 both Rq and 1/Q move 
towards each other, at a certain scale Q = 1/Rq they meet and then pass each other. 
Hence, the saturation occurs at higher values of than in a model with a fixed Rq. In 
addition the transition from high to low is faster. The line given by the condition 

^o(^) = ^ (16) 

will be called the critical line. In the parton picture it would describe the boundary of 
the critical density |[T^. The precise saturation pattern is determined by the fits of the 
three parameters in (§) to the existing DIS inclusive data. 



2.2 The Mellin Transformation 

In this section we briefiy describe the technique we use to perform the fits. In order to 
evaluate the cross section (0) we employ the Mellin transformation to factorize the wave 
function from the cross section. This reduces the number of numerical integrations to one 
and allows for an easier discussion of the scaling behaviour: 



^^'2 / ^ \ l+2iu 



0" X, r 



:i7) 



dv 



l/2+iv 



m. 



where 



m, 



X 



E 



^2 6ae^ £2(3/2 + zi.)r(l/2 + iv) ( + Am)^ 
27r ^ 



(1 + 3zz/)g2 + (3 + 2w)m] 



Q2 + 4m^ 
Q2 + Am] 



r(2 + w) 



Q2 + 



2i^i - + 



3 + 2iv 
4 



1 1 3 



2 



2' 2' Q2 + 4^2^ 



2 ' 2' 2' Q2 + 4^2^ 



and 



, \ _ 6aen^ 1^(3/2 + zi/) /g2^4m|\ 

^^^'''Q^j " r(2 + zz.) I Q2 ) 



/ 

16 p /3 ^ . 17 Q2 - 
""so ^^^2 '''2'2'g2 + 4^2^ 



(19) 



In addition in (|T^ 



G(^) 



rff2 (f2) 



-3/2-ij. 2n 
^('^ ) 



(20) 
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In order to have the right threshold behaviour and a smooth transition in the hmit 
— s> we also modify the Bjorken variable 

^ = ^ g2 • (21) 

The Mellin transformation for the functional form (|TUp of the function g which defines 
the dipole cross section in our analysis yields: 

G{y) = -n-^-zu). (22) 

For the sake of completeness we also present the corresponding expressions for the alter- 
native models in (|ri|) and (p!2D, 

^('^^ = T7^— fJr^ (23) 

1/2 + lu cosh(7rz/) 

and 

G(z^) = -T^ . (24) 

^ ' 1/2 -iu cosh(7ri^) ^ ' 

The main purpose of introducing a mass is to have a finite limit 0. However in 

the region Q ^ 1 GeV, where we mainly fit the data, the mass has little effect and one 
might as well consider the case with zero mass. In this case the functions (p!8D and (p!9[) 
reduce to 

( m _ \- 2 6ae™ vr 9/A + u' ( vr ^ sinh(7rz/) T{^/2 + iv) 
H^(i.,0) - ^Cf^- Losh(7rz/) J -T[-l/2-iv) 



and 



TT f n\ 2 Sftem TT 1/4 + z/^ / 71 ^ sinh(7r;/) r(3/2 + iz/) 

= ^ 8 TT^ l^^^ihM j -r(-i/2-w) ■ ^^^^ 



The results ( pSj) and ( pGj) are much simpler than in the massive case and probably familiar 
to those readers who have calculated the Mellin transform of the quark-box diagram in 
momentum space. The only difference is the ratio of the F functions at the end of each 
expression which is a residue of the Fourier transform. 

We use formulae (|1^ through (^2]) in our fits and perform the integration over u 



numerically. The main virtue of this representation is the quick convergence of the in- 
tegrand in ([T7|) which makes the numerical integration very fast. Another virtue is the 
rather simple analytical structure which allows to extract the leading scaling behaviour 
of (Jt for large as well as low Q^, as will be demonstrated in the next section. 
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2.3 Simplified Parametrization and Critical Line 



In this section we analyse the analytical structure of the cross section ( p!7D in the complex 
z/-plane. We use the formulae ( P^D and for the massless case. The position and 
characteristics of the singularities in the complex i/-plane determine the behaviour of 
the analysed cross section. In general the //-integration runs along the real axis, and 

depending on the argument ^ in eq. (|T7D one can close the contour in the upper or 

lower part of the complex plane. For example, in the case of large and not too small x 
("hard" regime) the mentioned argument is less than 1 and the contour has to be closed 
in the lower plane. The first singularity encountered is a pole at v = —i/2. Depending 
on the model for the dipole cross section given by (|2^), ( p3D or (p^), it can be a double or 
triple pole. The model (|2^ ) used in our analysis leads to a double pole which generates a 
logarithmic behaviour of the cross section flTi 



(27) 



One should note that the logarithm is due to the photon wave function (the quark-box 
diagram) and is related to the splitting of a gluon into a quark-antiquark pair. The factor 
in front of the logarithm arises because 1/2 + iz/ = 1 in eq. ( p!7[ ) and simply reflects the 
basic scaling behaviour of F2 combined with a certain power behaviour in x. One also 
recognizes that the longitudinal contribution ( P^ has only a single pole and therefore 
does not produce a logarithm. Since it is not the leading contribution we will ignore it in 
the following. 

In the "soft" regime where ( is larger than 1 we close the contour in the upper 

plane and find the leading pole at u = i/2. For the transverse cross section together with 
model (|22|) it is again a double pole leading to the following behaviour of flTi 



In 



X J 



{2i 



We can now combine the "hard" and the "soft" terms given by eqs. ( p7| ) and ( pH] ) into 



one expression 




+ 




In 



X 



Xr\ 



' Ql 



(29) 



where we have added 1 in the argument of the logarithms in order to allow for a smooth 
transition between the "hard" and the "soft" regimes. We have also introduced the 
parameters with a prime to indicate that the functional form should be refitted to get 
a good description. Although eq. (|29|) is a rather crude approximation of the original 
approach it reproduces the main features. 

As was discussed in our qualitative analysis we define the saturation as the transition 
from short distances to long distance with the characteristic scale given by the saturation 
radius Ro{x). Looking at eq. (p9|) we realize that the transition from the "hard" into 
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the "soft" regime occurs when ( ^ = 1. This equahty defines basically the same 
critical line as in our qualitative analysis of section |2.1| (see eq. ([16D ) 

(30) 




The precise location in the (x, (5^)-plane and the slope of the critical line is determined 
from the fits discussed in the following section. 



3 Discussion of the Fits 

We will discuss separately fits with and without charm contribution. As we will see charm 
has a quite strong influence on the fit and causes the critical line to move towards smaller 
scales. For the three light flavours we assume a common mass of 140 MeV which leads 
to a reasonable prediction in the photoproduction region. The general dependence of the 
total cross section a'^ ^ on the quark mass is such that the photoproduction cross section 
increases logarithmically with decreasing mass and diverges in the limit of zero mass. 
Thus the quark mass plays the role of a regulator for the photoproduction cross section. 

For our fit we use all available DIS-data below the Bjorken variable x = 0.01 combining 
systematic and statistical errors in quadrature. For the numerical evaluation we have 
implemented formulae (|l^ through (0). 



3.1 Light Flavours 

The result for our first fit with only light flavours is listed in the first row of Table 1. The 
corresponding plot in Fig. ^ shows only a subset of the data that were fitted. The most 
remarkable feature of the cross section in Fig. ^ is the turn over of the curves towards 
small values. This illustrates the change between the scaling and saturation region 
for high and low Q^, respectively, and is well reproduced by our model (solid lines). The 
critical line computed from eq. (^) is plotted across all curves and indicates the onset of 
saturation. For comparison we also show the corresponding cross section when the quark 
mass is set to zero (dotted lines). This line demonstrates that indeed the turn over occurs 
basically along the critical line and is not significantly influenced by the quark mass. 

Fig. ^ shows the logarithmic Q^-slope of F2 computed for fixed x and then plotted for 
fixed W"^ as a function of and x separately. These plots are inspired by an analysis 
which was carried out by ZEUS to stipulate the deviation from the conventional 
pQCD-approach at low values of (see also |^). The remarkable property of the 



presented plots is a distinct maximum for each of the slopes. The critical line lies slightly 
to the left of the maxima in both plots what might suggest that an alternative definition 
of the critical line could be introduced as being the path along the maxima. 

In Fig. 1^ we show the position of the critical line obtained in our analysis in the {x, Q^) 
plane. We observe that going along the critical line from x = 10~^ to a; = 10~^ we increase 
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FIT RESULTS 




(To (mb) A 


xo /{no.exp.po.) 


no charm 


23.03 


0.288 


3.04- 10-4 


441/372 = 1.18 


with charm 


29.12 


0.277 


0.41 • 10^4 


558/372 = 1.50 



the saturation scale from approximately 1 GeV^ up to 2 GeV^. This means that even at 
a future ep-coUider we expect only a rather small shift in the saturation scale. With great 
optimism it might go up to 3 GeV^ at a 1 TeV machine which, nevertheless, reaches quite 
far into the perturbative regime. 

An important prediction resulting from our analysis is the longitudinal structure func- 
tion Fl which is shown in Fig. |[ As we see Fl constitutes roughly 20% of F2 for 
around 10 GeV^ which is a reasonable value. 

We now concentrate on the W"^ dependence of the total cross section a'^ ^ for fixed 
Q^. We are particularly interested in the effective slope of a'^*^ as a function of W"^. This 
dependence is shown in Fig. ^ in a broad range of values and in Fig. | with a particular 
interest for smaller values of < 6.5 GeV^. Moving from high to low we see how the 
slope of a'^'*'^ flattens off indicating the low-Q^ saturation. One can also recognize a slight 
curvature in particular for = 0. This is associated with saturation in energy. 

A more explicit way of exposing saturation is achieved by plotting the effective slope 
Xeff computed from the relation 



d\n{W^) 



Kff = ^-^^ (31) 



as a function of for fixed values of x = (Q^ + 4mj) / {Q^ + W"^). This means that we 
move along the dashed lines shown in Fig. |^ when computing the effective slope. The 
resulting curve is shown in Fig. |^. The transition from the "soft" to the "hard" regimes 
is clearly visible. Quite interestingly, Ae// at very low turns out to be close to the 
standard value of 0.08 for the soft Pomeron. This value is dependent on the choice of 
the quark mass which is not completely constrained. However, if the mass is chosen 
such that the cross section in the photoproduction region roughly agrees with the data it 
automatically leads to the right value for the slope in this region. The other important 
feature of small-x saturation is that with decreasing x the slope Ag// also decreases, as 
can be inferred from Fig. ^ In the true high energy asymptotics it would approach zero 
which is related to the fact that in the complex plane of angular momentum j the leading 
singularity for our model is located at j = 1. 



3.2 Charm 

Charm gives a substantial contribution to the total cross section cr'*'*^ and cannot be 
ignored. We have performed a separate fit including charm but without introducing new 
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parameters. Technically this means that in the basic formulae (0-^ the sum over flavours 
has to be extended to include charm with a mass of 1.5 GeV. Also effected is x in the 
corresponding dipole cross section (]8|,pJ]) since it contains the quark mass (see eq. ([21|)). 
We do not change the light flavour mass and keep it equal to 140 MeV. The impact 
of charm is shown in the second row of Table. 1 and in the last three plots. The of 



the fit has increased but is still acceptable. This is also reflected in Fig. |T0| in which the 
agreement with the fitted data looks satisfactory. The most important change is a shift 
of the critical line down to smaller scales (see Fig. p^O]). 



More visible is the effect of introducing charm in the Q -slope of F2 shown in Fig. [11 
The shape of the peak becomes broader and the value of the slope at the high end of 
and X increases. The main reason for broadening is the fact that charm saturates at a 
much higher scale due to its high mass. We also see that the critical line moves further 
away from the maxima. We still believe that the critical line characterizes the dynamical 
part of saturation in contrast to the saturation for charm which is entirely due to its small 
size and predominantly perturbative. But even a small size contribution will at extremely 
large W"^ be bigger than the saturation radius Ro{x) and undergo small-a; saturation. 

A rough comparison with preliminary data on the charmed structure function shows a 
reasonable agreement with the data at low x. We have again looked at the effective slope 
Xeff, shown in Fig. |TT|, and find a slight increase at the low end of of about 10% as 
compared to our previous fit with light flavours only. 



4 Implication on Diffraction 

One of the important features of the wave function formalism is its flexible applicability 
for inclusive and diffractive scattering in DIS. We can immediately write down the cross 
section for diffractive scattering using the photon wave function defined in section El (see 



also ref. |18 



da,"" 



dt 



t=o 



^ JdhJ^'dz\^T,L{z,r)\'\aix,r')\' 



(32) 



and perform the Mellin transformation as in section 
= l^ldh£dz\^r,LMf 



daS 



dt 



dr 



12 



rl2 



\ l+2iv 



a(x,r'^)|^ (33) 



dv 
2^ 



X 



1/2+iv 



Ht,l z^, 



m. 



The only change appears to be the function G^{v). Provided the integrand in the second 
line of eq. (|33|) has at least a single pole at z/ = — i/2 we find the power behaviour with 
respect to the Bjorken variable x to be the same for the inclusive as well as diffractive 
cross section. The explicit calculation of the function defined by 



G^fz/) 



df' [f')-'"-'"' (1 - 
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= (2^/2+- _ 2) r(-l/2-iz/) (34) 
= (2 - 2^/2+*'^) G(z/) 

reveals that G'^(z/) is almost identical to G(i/) except for a factor which gives zero at 
V = —i/2. This zero has a common origin for all sensible models and simply reflects 
the fact that at small distances r the cross section vanishes proportional to which in 
diffraction turns into r^. In section p^ , eg. (pSj) , we found a double pole at z/ = 2/2 in the 



integrand for the transverse part which is now reduced to a single pole. Since we are left 
with at least a single pole, the diffractive 7p-cross section for transverse polarized photons 
has indeed the same power behaviour in x as the inclusive 7p-cross section. 

We also found in section pl2| , eg . (p6D , that the inclusive longitudinal cross section 
had only a single pole which disappears completely in the diffractive cross section. The 
conclusion is that the next pole at z/ = —3i/2 becomes the leading pole and the power 
1/2 + 11^ in eg. (|33| ) eguals 2. Conseguently, the longitudinal part in diffractive scattering 
rises at small x two times as fast as in the inclusive case. Our result for the longitudinal 
contribution also confirms the well known higher twist nature of the diffractive cross 
section. 

It is important to note that the results discussed in this section are derived without 
imposing cuts on the final state except for the basic rapidity cut. In particular, there has 
no cutoff being assumed for the transverse momentum of the final state partons. This 
allows to use formula (P^). In ref. p2[ a similar study lead to a closely related statement 



about the conservation of "anomalous dimensions" . In an extended version with non-zero 



momentum transfer links can be found to conformal invariance in the Regge limit |^ . 

We can go a step further and calculate the ratio of the diffractive and inclusive cross 
section in the limit of large Q^. This can be done by shifting the contour of the u- 



integration down the lower half of the complex plane, as was done in section |2.3| , and 
picking up the leading pole: 

^ ^ ^0 3 ln(2) 

atot 4:7rB 6 In + 6 A ln(a;/xo) + 6 + 1 ' 

where B denotes the diffractive slope. A mild logarithmic dependence on x is found in 
the ratio which gives a slight logarithmic rise with decreasing x. Increasing we find 
a logarithmic suppression which is due to the extra zero in G^(z/) aX u = —i/2. Taking 
the parameters found in our fit and the experimental value for B we find that the ratio 
turns out to be egual to a few percent. We have to note that a contribution of similar 
size is expected from the emission of a gluon. This contribution will have the same power 
behaviour in x as for the guarks and we expect a ratio of a similar size. All contributions 
together will be of the order of 10%. 

More precise studies on diffraction based on the model discussed here will be presented 
elsewhere p4 |. 
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5 Summary 



We have presented a model which provides a good description of all DIS data below 
X = 0.01 (including the photoproduction data). An important ingredient of our approach 
is the presence of small-x saturation. This is achieved by introducing an x-dependent sat- 
uration radius Ro{x) = (x/xq)^^"^ which we define with the help of two parameters. These 
parameters, together with the overall normalization of the cross section, was determined 
by fits to the DIS data. The parameterization obtained in this way was then extended 
into the photoproduction region showing reasonable agreement with the data there. A 
particular result of this extrapolation is the effective Pomeron intercept of approximately 
1.08. In the DIS region this effective intercept increases to the asymptotical value of 
1.29. We have introduced a "critical line" which marks the boundary of the saturation 
region in the {x, Q^)-plane. Finally, we found an intriguing result when applying our 
model to diffractive scattering processes, namely, an almost constant ratio of diffractive 
and inclusive DIS cross section. 

Our model is basically phenomenological. It is remarkable, however, that with only 
three fitted parameters we successfully described the small-x data. The theoretical basis 
is the separation of the perturbative dissociation of the photon into a quark-antiquark 
pair and the subsequent interaction of the pair with the proton. The latter is modelled 
in an eikonal way. The main consequence is that we do not only have saturation at 
low Q^, but also saturation at low x. The natural extension of our model would be the 
incorporation of a perturbative treatment for large Q^. Also important is the development 



of a microscopical picture behind saturation as proposed in ref. p5[ . 

A real test for our model would by a future ep-coUider. Our prediction is that with 
increasing energy the saturation scale moves up in Q^. For an 1 TeV collider it woud be 
roughly a factor of 2. 
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Figure 3: The 7*p-cross section for various energies. The sohd hues show the fit results 
with a light quark mass of 140MeV. The dotted lines show the cross section with the same 
parameters but with zero quark mass. The line across the curves indicates the position 
of the critical line. 
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Figure 4: The logarithmic Q^-slope of F2 for fixed energies W, plotted as a function of 
Q'^ and x. The line across the curves shows the position of the critical line. 
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Figure 5: The position of the critical hne in the [x, (5^)-plane. The narrow hatched area 
corresponds to the acceptance region of HERA. The wide hatched region indicates the 
range for a future 1 TeV ep-colhder. The boundaries are hues of constant y. 
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Figure 8: The same cross section as in Fig. |^ with the emphasis on small-Q^ values. 
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Figure 9: The effective slope in W"^ of the total cross section, Ag//, plotted as a function of 
for fixed x. The solid line appears twofold due to two different methods of calculation. 
One is based on the analytical expression (0), the second is calculated from fig.0 by 
taking a pair of points along the lines of constant x = 10~^ and 10~^. 
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Figure 10: The 7*]9-cross section including charm. The charm contribution itself is plotted 
as dotted line. 
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Figure 12: The effective slope Ae// based on the charm fit. The intercept is shghtly higher 
at low Q"^ as compared to fig.pl. 
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